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SHOALING  OF  NONLINEAR  WAVE-GROUPS  ON  WATER 
OF  SLOWLY  VARYING  DEPTH 


Ruth  Iusim  and  Michael  Stiassnie 

Department  of  Civil  Engineering,  Technion  I.I.T. , 
Haifa  32000,  Israel 


Abstract 

An  approximate  analytical  solution  describing  the  shoaling  of 
modulated  wave-trains  is  presented.  This  solution  provides 
new  information  about  the  wave  field  evolution  as  well  as 
about  the  wave— induced  mean  current  and  set  down. 


1.  INTRODUCTION 

The  shoaling  of  weakly  nonlinear  surface  wave  groups  is  important 
to  the  understanding  of  coastal  wave  climate  and  coastal  flow  regime. 

In  the  past,  most  efforts  concentrated  on  the  equally  important 
though  Ampler  problem  of  shoaling  of  wave-trains  (i.e.  monochromatic 
wave  groups),  for  details  see  Stiassnie  &  Peregrine  (1980). 

The  first  mathematical  formulation  for  shoaling  of  wave-groups 
was  given  by  Djordjevic' and  Redekopp  (1978),  and  in  a  somewhat  improved 
version  by  Stiassnie  (1983).  This  formulation  is  limited  to  cases 
where  the  water  depth  is  small  compared  to  the  group-length.  Equations 
suitable  for  water  depths  of  the  order  of  the  group-length  are 
deduced  in  Peregrine  (1983);  combining  the  constant  depth  model  by 
Davey  and  Stewartson  (1974)  and  the  higher-order  model  for  infinitely 
deep  water  by  Dysthe  (1979). 

The  only  available  solutions  are  those  for  the  shoaling  of  isolated 
wave-packets  (solitons) ,  which  were  originally  given  by  Djordjevic' 
and  Redekopp  in  their  1978  paper.  They  predict  that  a  soliton 
envelope  can  undergo  fission  only  if  it  propagates  into  deeper  water. 

By  heuristic  assumptions  for  the  evolution  along  the  slope,  they  also 
estimate  the  number  of  solitons  emitted  after  a  single  soliton  descends 
from  a  shallower  shelf.  A  more  recent  study,  Turpin,  Benmoussa  and  Mei 
(1983)  confirms  these  results  qualitatively,  but  not  quantitatively. 


To  the  best  of  our  knowledge,  no  results  for  shoaling  of  wave- 
groups  (i.e.  modulated  wave-trains)  have  been  presented  so  far.  These 
modulated  wave-trains  are  of  particular  importance  since  almost  every 
wave-train  will  eventually  become  modulated  due  to  its  intrinsic 
Benjamin-Feir  instability.  The  aim  of  the  present  paper  is  to  throw 
some  light  on  the  evolution  during  the  shoaling  of  a  modulated 
wave-train  and  its  influences  on  the  mean  free  surface  and  the  wave- 
induced  mean  flow. 

Sections  2,  3  and  4  outline  the  derivation  of  the  mathematical 
model  and  its  simplifications  to  a  level  which  enables  an  analytical 
solution.  The  model  presented  at  the  end  of  section  2  is  rather 
general;  it  consists  of  a  coupled  system  of  equations  for  the 
complex  wave  envelope  and  the  induced  mean  flow  potential.  In  section 
3  we  add  the  assumption  of  periodicity  of  the  modulation,  which  leads  to 
decoupling  of  the  system  of  equations.  The  result  is  a  nonlinear 
Schrddinger  equation  with  a  variable  (depth  dependent)  coefficient, 
which  can  be  solved  by  reasonable  numerical  efforts.  In  section  4 
we  adopt  the  assumption  of  three-waves  systems,  which  has  been  used 
for  constant  depth  in  Stiassnie  and  Kroszynski  (1982),  and  which 
together  with  a  W.K.B.  type  approach  enable  an  analytic,  though 
asymptotic  solution.  This  asymptotic  solution  is  compared  with 
numerical  results  in  section  5.  Sections  6,  7  and  8  include  a 
detailed  presentation  and  discussion  of  the  physical  results  obtained 
from  out  calculations. 


2.  EVOLUTION  EQUATIONS 


Assuming  irrotational  motion,  there  exists  a  velocity  potential 
$(x,z,t)  which  satisfies  Laplace's  equation: 


<j>  +  <(> 

xx 


zz 


0, 


(2,1) 


t  is  time,  x  is  the  horizontal  coordinate  in  the  direction  of  wave 
propatation,  and  z  is  the  vertical  coordinate  pointing  upward  from 
the  undisturbed  free  surface. 


The  boundary  condition  on  the  bottom,  z  =  -h(x),  is 

♦z  “  -h' <x)*x.  (2.2) 

The  boundary  conditions  on  the  free  surface,  z  =*  C(x,t),  are  the 
kinematic  condition: 

♦.  -  't  +  ♦r'x*  <2‘3) 

and  the  dynamic  condition: 

2gC  +  2 A  4-  +  <P^  *  0-  (2.4) 

Considering  situations  for  which  the  depth  h(x)  as  well  as 
the  wave  input  vary  slowly,  we  assume  all  wave-field  properties 
to  be  slowly  changing.  To  render  the  term  'slowly'  explicit  we 
introduce  a  small  nondimensional  parameter  c  which  is  a  measure  of 
the  wavy  surface  slope,  and  define  the  following  new  variables: 


-4- 


where  ft*  *  3ft/3K  is  the  group  velocity.  ft2  =  gKth(Kh)  is  the  linear 
dispersion  relation,  relating  the  leading  order  constant  frequency 
ft  to  the  leading  order  wave-number  K(^)  and  the  water  depth  h(^). 


Having  the  wave-groups  (i.e.  wave  fields  with  narrow  spectra) 
in  mind  the  velocity  potential  cf>  and  the  free  surface  displacement  £ 
are  expanded  in  Fourier  series: 


■f  0  2,*i  0 

<(>  *  $o(x,z,t)  +  +  4>2(x,5,z)e  +  ...  c.c.}  (2.6a) 


5  =  50(t,5)  +  {C1(t,5)«10  +  C2(T,5)e219  +  ...  +  c.c  >  (2.6b) 


a 

where  9  *  (/  K(£)dx  -  ftt) ,  and  c.c.  stands  for  the  complex  conjugate. 
With  e  chosen  to  be  small,  the  functions  <^(r,£,z)  and  ( t,£>  for 

j  >  1  are  expanded  formally  in  power  series  of  e  : 


‘J’jCr.s.z)  «  E  en*  (x,g,  z) 
J  m*j 


(2.7a) 


“  2  e,tt4mj(T»£> 
m*j  J 


(2.7b) 


The  mean  water  level  and  the  induced  mean  flow  potential  <j>o> 
require  a  special  treatment  and  are  best  written  as: 


«0(t,5)  -  e2C2Q.  (2.7c) 

4»o(x,z,t)  *  e$10(t,£,z)  +  E2<^20(5)*t.  (2.7d) 

and  £«j 


Note  that  it  turns  out  that  while  <£ 


for  j  >  1  are  all  of 


> 
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order  one,  the  order  of  ^Qf^o  ^20  *"s  ^etween  one  f°r  ^  =  0(1) 
and  e  for  Kh  -*■  00  .  The  term  s2  ^^Q(C)*t,  in  Eq.  (2.7c),  is  needed 
to  suppress  terms  that  grow  boundlessly  with  time  at  higher  order, 
as  was  shown  in  Stiassnie  (1983).  Following  the  method  of  derivation 
used  by  Djordjevic'  and  Redekopp  (1978)  but  using  Eqs.  (2.7c,d), 
instead  of  expanding  <j>Q  and  5  in  power  series  of  e  (which  is  justified 
only  for  Kh  =*  0(1)),  we  obtain  a  system  of  evolution  equation  (2.9a,b,c,d) 
for  and  the  complex  wave  envelope. 


as  follows: 

i  aa*  . 


A(t,5)  = 

-2i^ll 

(2.8) 
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(2.9a) 
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-h  ^  z  <  0, 

(2.9b) 

^lo  +  g“  *10 

Z  °  TT 

g&2  ,  ,0 

~20-  <IAI  >t 

;  z  =  0 

(2.9c) 

♦10,  * 

z=-h. 

(2.9d) 

The  depth  dependent  coefficients  are: 

,  9- 12th2 (Kh)+13th4 (Kh)-2th6 (Kh) 

S1  “  20  *  8th3 (Kh) 


S2  ■  H  ■  &  +  sech2  (Kh) ) 


03  -  2jj-  sech2  (Kh) 


(2.10a) 

(2.10b) 

(2.10c) 


A  set  of  modulation  equations  equivalent  to  (2.9),  but  for  constant 
depth,  has  been  recently  derived  from  the  finite  depth  Zakharov  equation 


see  Stiassnie  and  Shemer  (1984). 


The  mean  water  level  C 

o 


is  given  by: 


sech2(Kh)  •  |  A|  2 


(2.11) 


and  is  of  order  between  e2  to  c3  ,  depending  on  the  water  depth. 


3.  PERIODIC  CASES 


Restricting  the  discussion  to  cases  for  which  the  complex  wave  envelope 
A(t ,£)  is  periodic  in  t,  and  assuming  zero  averaged  (over  t)  mass 
flow  in  the  x-direction  (as  in  the  case  of  an  impermeable  beach)  enables 
the  decoupling  of  the  system  (2.9).  For  these  cases  A  is  governed  by 
the  nonlinear  SchrHdinger  equation: 


i 

2fl' 


an* 


A  +  iA^  + 


nr 


e^a. 


2(flT) 


A  - 

TT 


n' 


|a|2a  = 


lAp  A,  (3.1) 


where  the  bar  denoted  averaging  over  t,  and 

82B2 

2n(i-gh/(n’)2) 


“l  “  si  + 


(3.2a) 


s282 

ZBL  (a  +8  )  .  s 

2hfi  2a(l-gh/(a,)z) 


(3.2b) 


The  induced  mean  flow  potential  is  given  by: 


where  Y  is  the  period  of  A(t£); 


a  =  cosh(2iren(z+h)/&,y) 


n  2Q{^  Sinh(^£ii)-C^ coshC^h) }  ~n> 
eQ’y  Qry  y  Q'y 


and  C  are  the  Fourier  coefficients  of  (IaI2-)  : 

n.  ’  '  t 


(|a|2)  -  S  tC  (rle2*lnT/T«  (5>e-2’lnT/Y> 

T  ,  EL  *  — n 

n=l 


(3.4) 


The  convergence  of  the  Fourier  series  (3.4)  was  assumed  to  be  inde¬ 
pendent  of  e  . 

The  potential  which  is  needed  to  calculate  the  mean  water  level 
£q,  see  Eq.  (2. 11),  is  given  by: 


(f>2Q  — e  “|A|Z  gKftV2flh. 


(3.5) 


Hq.  (3.1)  is  identical  to  Eq.  (4.9)  in  Stiassnie  (1983),  which  has 
been  derived  assuming  no  vertical  dependence  of  the  Induced  mean  flow. 
The  latter  assumption  is  strictly  valid  only  when  Kh  =  0(1),  namely, 
for  cases  where  the  water  depth  h  is  small  compared  to  the  group  length 
(er)  *.  For  these  cases  the  present  $  and  £  yield  Eqs.  (3.7) 
and  Eq.  (3.8)  of  Stiassnie  (1983). 


A  simpler  dimensionless  form  of  Eq.  (3.1)  is  obtained  by  means  of 


the  transformations 
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-l#2QSa\3i  .  ,  ,Tl rrr  rX  a2dXs  . 

“  e  ( - 3”)  A  i  (lAl  J  ~rp~^  » 


T  =  x/y; 
which  give 


f 


Y‘  '  2(0') 


TTTdC. 


(3.6a) 

(3.6b) 


i'J'v  +  <J»TT  +  u[*l2+  =  0, 

1 


rX  '  TT 

-g3a'y2a 

U(f)  - 


5o  " 


0=0 


(3.7) 

(3.8) 


The  dimensionless  parameter  u  is  a  monotonic  increasing  function 
of  Kh,  having  the  values  zero  and  (Qy)2  for  Kh  =  1.363  and  Kh  -*■  ® 
respectively.  The  statement  of  the  mathematical  problem,  given  by 
Eq.  (3.7),  is  completed  by  the  following  input  condition  at  X  *  0 
(i.e.  x  *  x  *•  a  reference  point  in  infinitely  deep  water). 


iKT.0)  =  1  +  28ei0tcos(2irT) 


(3.9) 


which  corresponds  to  a  system  composed  of  a  carrier-wave  and  a  symmetric 
"side-band"  disturbance 


Re{e-lflt+ee-i|:a+2''e/nT)!lt-*d+ 


+  Be-lCa-2.e/aT>!!t-aJ)+<|(c2) 


(3.10) 


For  constant  depth  ,  u  *  const,  it  is  well-known  that  Eq.  (3.7)  with  Tin  (0,1), 
subject  to  periodic  boundary  conditions  has  the  following  X  invariants 
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M 

MM 


•  V 


w 


mm 


•-  a 
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(3.11a) 


Jj  -  /  m  dT 


(3.11b) 

(3.11c) 

These  invariants  are  determined  by  the  inout  condition,  Eq.  (3.9)  so 
that 

J1  =  1  +  262  ;  J2  -  0  (3. 12a, b) 

J3  =  1  +  (4-P+2cos2a) • 282+  664  (3.12c) 

where 

P  «  8tt2/u  (3.13) 

For  varying  depth,  u  =  u(X),  and  J2  remain  invariant  and  are  given 
by  C3.11a,b)  and  (3.12a,b),  but  is  a  function  of  X  governed  by  the 


J2  =  I  0WT  -W»*)dT 
o 


1 


J3  -  /  (kl4  -  kT[2)dT 
o 


equation. 
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4.  THREE- WAVES  SYSTEMS 


The  solution  of  Eq.  (3.7)  can  be  expanded  in  a  Fourier  series 


<KT,X)  =»  E  D  (X)e2TCinT  (4.1) 

n— “  n 

The  boundary  condition  at  X  =  0,  Eq.  (3.9),  gives  Dq(0)  «*  1; 

Dj  (0)  =*  D_1(0)  -  Be1*;  \(0)  -  D_fl(0)  =  0  forn  >,  2. 

Stiassnie  and  Kroszynski  (1982)  truncated  the  above  given  series  and 
considered  only  three  waves  systems: 


iKT.X)  -  y.  D  (X)e27rinT 

i  a 

n=-l 

Substituting  Eq.  (4.2)  into  Eq.  (3.7)  yields 
ordinary  differential  equations: 


(4.2) 

the  following  system,  of 


dD 

idr  +  ^(lD0l2+4iDii2)v2Di^  “  °»  (4.3a) 


dD,  p 

+  uC(2[Do|2+3|D1|2-  |)Dj  +  Dj)*]  -  0.  (4.3b) 

Note  that  Eq.  (3.11b)  yields  D_^  =  D^. 

For  constant  depth  the  system  of  Eqs.  (4.3a,b)  has  exact  solutions 
in  terms  of  Jacobian  elliptic  functions  with  periods  of  order  1  in  X  which 
is  summarized  in  the  Appendix;  for  details  see  Stiassnie  and  Kroszynski 
(1982).  These  solutions  depend  on  the  invariants  and  on  the  parameter 

U  ,  which  in  turn  depends  on  the  water  depth  h  and  on  the  modulation 
period  y  .  For  very  mild  depth  variations,  where  h„  *  o(l),  we  apply 

A. 


an  asymptotic,  WKB  related  approach,  assuming  the  local  solution  to  be 
that  of  the  constant  depth  type  and  using  Eq.  (3.14)  to  determine  J^. 

and  y  are  fixed  by  the  input  conditions  and  J^,  is  given  through 
I(P)  by: 


where 

I(P)  =  J3(P)  -  (4.5] 

The  initial  value  of  I,  at  X  =>  0,  where  P  »  P  is  denoted  by  I  and 

o  o 

is  given  by: 

I  -  &2[282  +  4(l+cos2o)-2Po]  (4.6) 

Iq,  as  well  as  I(P)  were  assumed  to  be  of  o(l)  throughout  the  rather 
lengthy  derivation  of  Eq.  (4.4).  In  all  our  examples  we.choose(Y  =*  2irft 
PQ  *  2,  corresponding  to  the  fastest  growth-rate  of  the  Benjamin-Feir 


instability. 


5.  NUMERICAL  VERIFICATION  OF  THE  ASYMPTOTIC  SOLUTION 

In  order  to  appraise  the  relevance  of  the  asymptotic  solution  given 
in  the  previous  section,  we  compare  its  results  with  those  of  a  numerical 
solution  of  the  system  of  ordinary  differential  equation  (4.3a,b). 

Fig.  5.1  shows  I=I(P)  for  four  initial  values  of ^=1(2) (IQ  =  -0.04, 
-0.01,  0.04,  and  0.1).  The  broken  line  represents  the  asymptotic  solution 
and  was  obtained  by  numerical  integration  of  Eq.  (4.4).  The  solid  line  was 
obtained,  by  substitution  of  the  results  obtained  from  a  numerical 
solution  of  the  system  of  O.D.E  (4.3a.b)  into  the  expression 

I(P)  =  2|D1i2{|D1|2+2|Do|2-P+2|DQ|2cos[2(arg  D^arg  Dq)]}  (5.1) 

Fig.  5.1  about  here 

The  numerical  solution  of  the  system  of  O.D.E  was  obtained  using 
a  trapezoidal  method  and  assumign  the  P(h(X))  =  2+0. 2X.  Note  that  the 
assumption  |l|  «1,  which  is  necessary  for  the  asymptotic  solution  to  be 
valid,  imposes  a  restriction  on  the  range  of  variation  of  P  (P  =  2.8 
corresponds  to  02h/g  ~  4).  In  Fig.  5.2  we  show  three  parts  of  the 
exterior  group  envelope  |$(0,X)|  as  well  as  the  interior  group  envelope 
|*(Js,X)|  for  the  input  conditions  cc  »  0,  8  »  0.158  (IQ  ■*  0.1). 

Fig.  5.2  about  here 


Here  again,  solid  lines  represent  the  numerical  solution  of  Eqs.  (4.3a,b) 
with  P  *  2+0. 2X  while  the  broken  lines  correspond  to  results  obtained 
by  the  asymptotic  method,  utilizing  the  relation: 


[4>(T,X)[2  »  {-21+2 (4J L-P) z-7z2+  [s(2I+2Pz-z2)'i+4zcos2irT]2  }/8z  (5.2 


where  z  is  given  in  the  Appendix. 

The  three  parts  shown  in  Fig.  5.2  are  for  P=2,  2.44,  2.75  for  the 
asymptotic  solution  compared  to E  in(2.,2.Q3),  (2.44,  2.49),  (2.75,  281) 
for  the  numerical  solution  of  the  O.D. E. , respectively.  The  agreement 
between  the  two  methods  of  solution,  as  seen  in  both  the  above  figures 
is  rather  encouraging  and  seems  to  indicate  the  validity  of  our  new 
asymptotic  solution  of  the  system  (4.3a,b).  Nevertheless,  one 
still  has  to  answer  the  question  if,  and  to  what  extent,  the  syetem 
(4.3a,b)  itself  is  a  reasonable  substitute  for  the  N.L.S,  Eq.  (3.7) 
For  constant  depth  Stiassnie  and  Kroszynski  (1982)  show  a  good 
quantitative  agreement  in  the  length  of  the  modulation-demodulation 
cycle  and  only  a  qualitative  agreement  for  the  amplitudes.  A  similar 
trend  can  be  seen  in  Fig.  5.3,  which  compares  two  numerical 
solutions,  one  for  the  N.L.S.  (3.7)  -  dotted  line,  and  the  other 
for  the  system  of  O.D.E.  (4.3a,b)  -  solid  line. 

Fig.  5.3  about  here 

The  input  data  in  Fig.  5.3  is  a  *  0,  8  *  0.1  (Iq  *  0.04)  and  the 

variation  P  *  2  +  0.2X  is  assumed;  Both  the  exterior  and  interior 
group  envelopes  are  drawn. 


1 


We  believe  that  our  much- simplified  asymptotic  solution  is  not  over¬ 
simplified,  and  is  able  to  produce  quite  a  few  results  of  qualitative, 
and  maybe  even  serai-'viantitative  relevance,  which  enable  us  some  new 
physical  insight. 

6.  ON  P  AND  I  - 
o 

One  fundamental  property  of  the  asymptotic  solution  is  that  it 

deoends  on  P  and  I  solely.  Given  the  input  data  a,B(and  P  =2),  I 
o  '  o  o 

is  determined  by  Eq.  (4.6).  Then,  integrating  Eq.  (4.4)  from  Pq  to  P 
the  parameter  I(P)  is  found,  and  the  solution  given  by  Stiassnie 
and  Kroszynski  (1982),  is  locally  applied. 

Fig.  6.1  about  here 

Fig.  6.1  gives  the  relation  between  P  and  nondimensional  local 

water  depth.  K  h,  (where  K  *  ft2/g  is  the  wave-number  in  infinitely  deep 

wate$,  for  y  *  2rft  ?  is  a  monotonic  decreasing  function  having  the 

values  infinity  at  K  h  =  1. 195(Kh=l.  363)  and  2  for  Kh  +  ®  .  The  input 

value  Iq,  (for  Pq=2)  dependence  on  a  and  0  is  shown  in  Fig.  6.2. 

Note  that  different  combinations  of  a  and  0  give  the  same  I  ,  and  thus 

o 

basicly  the  same  solution  for  any  P. 


Fig.  6.2  about  here 
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7.  GROUP  ENVELOPES 


The  free-surface  of  an  (unstable)  shoaling  wave-train,  displays 
three  distinct  length  scales:  X^  -  the  wave  length;  X ^  -  the  modulation 
or  group  length;  and  X^  -  the  modulation-demodulation,  group-envelope, 
or  maybe  best  Supergroup1  length.  These  three  lengths  are  given  by 


XL  -  2ir/K, 

X^  =  2Tr£-14V/fl  , 

.  ,  4e''2P(fl')3  2P2(4-P)2 

3  Q2Q"/P(4-P)  (2P-1)I 

see  Fig. (7. 2a) . 


(7.1a) 

(7.1b) 

(7.1c) 


It  can  easily  be  seen  that  in  the  range  of  depths  where  the 
asymptotic  solution  applies,  K^h  p  4,  X^  and  X^  remain  almost 
constant.  On  the  other  hand,  X^ ,  which  depends  on  P,  exhibits  quite 
a  remarkable  variation,  as  shown  in  Fig.  7.1. 


7. 1  about  here 


In  Fig.  7.1  we  show  the  variation  of  X^  as  a  function  of  the  depth 
K^h  for  four  different  input  data  I  *  -0.04,  -0.01,  0.04,  0.1.  For 
IQ  <0,  Xj  decreases  with  decreasing  depth,  but  for  Iq  >  0  X^  increases 
with  decreasing  depth  up  to  a  "critical  depth"  (corresponding  to 
1*0)  and  from  there  on  starts  to  decrease. 


Fig.  7.2  shows  the  group  envelopes  (dashed  line)  and  wave 
envelope  (solid  lines)  at  a  fixed  instant  for  e  *  0.2,  at  the  following 


four  locations:  (a)  -  infinitely  deep  water,  P  =2,  I  =0.1; 

o  o 

(b)  K^h  =  11.2,  P  =  2.2,  I  =  0.052;  (c)  Kh  =  5.7,  P  =  2.45,  1  =  0.01; 
(d)  Kh  =  4.2,  P  -  2.75,  I  =  -0.028. 

Fig.  7.2  about  here 

In  Fig.  7.2a  we  have  added  a  portion  of  the  wavy-surface  (thin  solid 
line)  as  well  as  the  lengths  ^3*  Note  that  the  supergroups 

(namely:  the  exterior  and  interior  group  envelopes)  are  fixed  in 
space,  while  the  wave  envelope  moves  with  group-velocity  and  the  waves 
themselves  with  the  phase  velocity.  Similar  sketches  to  Fig.  7.2 
were  obtained  for  the  three  other  cases  given  in  Fig.  7.1. 

In  order  to  complete  the  picture  for  shallower  water  depth  we 
present  in  Fig.  7.3  the  group  envelopes  at  five  locations:  P  =a:(2,2.05) 
b: (2. 32,  2.63),  c:(17.7,  22.8),  d: (-22. 4,-17. 5) ,  e: (-4. 1,-3. 9) , as 
obtained  from  a  numerical  solution  of  the  system  (4.3a,b)  for  the 
same  input  data  as  in  Fig.  7.2  assuming  u  *  39.47-10X. 

Fig.  7.3  about  here 

The  results  in  Fig.  7.3  indicate  that  continues  to  shorten 
and  that  the  intensity  of  modulation  decreases. 
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8.  THE  MEM  FLOW  FIELD 


We  express  the  mean  flow  u  =3<j>  /3x  ■  -3tJ>/3z,  v  ~3<f>  /3z=3i|>/3x 

o  o 

through  the  stream  function  t|; 

’  =  (|doI^Idx|^  + 


g3Y C~L 


g0(D*D.-+D  Df ) sh(2 (Z+H) ) cos  (2ttT) 

Z  O  1  O  i 


tofr  s!i(2H)-ch(2H) 


?  -1  BJd.  [£h(4(z+H))cos(4irT) 

£  T e  ^  A 


4xflxft' 


sh(4H)-2ch(4H) 


+  constant  (8.1) 


—  - 

where  Z  “irez/ft'y  ,  H  ^irEh/O'y  and  “e  il3!(i/g2  are  dimensionless 
quantities.  The  constant  in  Eq.  (8.1)  is  chosen  so  that  f  =  0 
at  the  bottom.  The  mean  free  surface  is  given  by  Eq.  (2.11) 


Figures  3.1,  8.2,  8.3  about  here 


The  stream-function  'KTjZ)  as  well  as  the  mean  free-surface 
for  cases  a,  b  and  d  of  Fig.  7.2  are  presented  in  Figures  8.1,  8.2 
and  8.3  respectively.  These  figures  demonstrate  the  rather  complicated 
structure  of  the  wave  induced  mean  flow  field. 


Some  of  the  main  features  are:  (i)  the  mean  current,  which  is 
shown  in  part  (b)  of  the  figures,  a a  well  as  the  mean  free  surface, 
in  part  (c)  exhibit  a  somewhat  cellular  structure  influenced  by  the 
wave  envelope  variations;  (ii)  a  dominant  adverse  cu  rent  appears 
underneath  the  high  waves  and  a  much  weaker,  positive  current  (in  the 
wave  propagation  direction)  under  the  low  waves,  for  the  shallower  cases 
the  positive  currents  almost  disappears;  (iii)  the  magnitude  of  the 
maxi mum  adverse  currents  at  the  free  surface  is  almost  the  same  for 


all  three  depths  (KJi  =<*>,  11.2  and  4.2);  Civ)  one  can  notice  the 
tendency  of  the  flow  fields  to  become  more  uniform  in  their  lower 
parts  and  on  the  sides  of  the  supergroups  (where  the  modulation 
amplitudes  get  much,  smaller) ;  (v)  there  is  a  set-down  in  the  mean 
free  surface  accompanying  the  peaks  of  the  wave  envelope  and  a  smaller 
set— up  accompanying  their  troughs. 


APPENDIX  -  Solution  of  Eqs.  (4.3a,b)  for  constant  depth 


Input  data:  JpJj.P;  I  =  J  -J^ 


ci  =  yCAJ^p)  {l-  [l- 

Ul  "fa 

(A.  1) 

(AJj-P^  J  1 

c2  -  p  *0(1  +$)hl 

(A.  2) 

c  *  maxCCj^.C^) ; 

d  =  min(C1,C2) 

(A. 3a, b) 

a  =*  2P; 

b  «  yCA-P) 

(A. 3c, d) 

»  1  _  a~k  .  c~d 

a  b 

(A.  4) 

y.-  .  X 

(A.  5) 

cd(y,k)  is  a  Jacobian  elliptic  function  of  the  argument  y  with  modulus  k 


-  a»b« (1-cd2) 
a-b  *  cd2- 


(A.  6) 


fa[  -  1^*72;  |D0|  -  /j^i 


(A. 7a, b) 


cos[2  (argD  ^-argDQ)  ]  =• 


I-l.Sz^+CP-ZJjJ'z 
2z (J.-z) 


(A.7c,d) 
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Numerical  solution  of 


5.1  I=I(P),  -  Asymptotic  solution, 

Eqs.  (4.3a,b). 

5.2  Exterior-group-envelope  jijj(0,X)|  and  Interior-group  envelope  U(?s,X)  | 

for  the  input  conditions  a  *  0,  8  =  0.158,  (1^=0. 1)  -  Asymptotic 

solution,  -  Numerical  solution  of  Eqs.  (4.3a,b). 

5.3  Exterior-group-envelope  |ijj(0,X)|  and  Interior-group-envelopej^j^x)  |  for 

the  input  conditions  ci=0,8  =  0.1,  (1^=0. 04),  - Numerical 

solution  of  O.D.E  (4.3a,b)  Numerical  solution  of  N.L.S. 

(3.7). 

6.1  P=P(K  h)  for  P  =2. 

oo  O 

6.2  I  =1  (a, 8)  for  P  =2  (on  S2=  0:  I  « 0 )). 

o  o  o  o 

7.1  The  'supergroup'  length  X^  as  a  function  of  depth  h. 

7.2  The  group  envelope  ( - )  and  wave  envelope  ( - )  at 

t  -  constant  for  Pq=2,  Io=0. 1, G  =0.2  at  K^h  =  (a)  ®,  (b)11.2, 

(c)5.7  and  (d)4.2. 

7.3  The  group-envelopes,  for  P  =2,  I  =0.1,  at  K  h  -  (a)  (b)  7.11-4.8, 

o  o  00 

(c)1.32-1.29,  (d)  1.12-1.11,  (e)  0.96-0.95. 

8.1  The  flow  field  for  P  =2,  I  =0.1,  e  =0.2  at  K  h  *«  .  (a)  group- 

o  o  00 

envelopes  and  wave  envelope  (b)  mean  flow  stream-lines,  (c)  mean 
free-surface. 

8.2  As  in  8.1  for  K  h  =  11.2 

OO 

8.3  As  in  8.1  for  K  h  =  4.2 
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